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The microscopic structure of the prismatic S3 (1010) twin boundary in 
a-Al203 is characterized theoretically by ab-initio local-density- functional 
theory, and experimentally by spatial-resolution electron energy-loss spec- 
troscopy in a scanning transmission electron microscope (STEM), measur- 
ing energy-loss near-edge structures (ELNES) of the oxygen i^- ionization 
edge. Theoretically, two distinct microscopic variants for this twin inter- 
face with low interface energies are derived and analysed. Experimentally, it 
is demonstrated that the spatial and energetical resolutions of present high- 
performance STEM instruments are insufficient to discriminate the subtle dif- 
ferences of the two proposed interface variants. It is predicted that for the cur- 
rently developed next generation of analytical electron microscopes the pris- 
matic twin interface will provide a promising benchmark case to demonstrate 
the achievement of ELNES with spatial resolution of individual atom columns. 



I. INTRODUCTION 

The local atomistic and electronic structures at extended defects in materials, such as 
grain boundaries in polycrystals or interfaces in layered nanostructures, are highly influential 
on the structural or functional properties of many technological devices, and hitherto, strong 
efforts have been made to characterize them microscopically. Recently, significant progress 
has been achieved by combining experimental imaging and spectroscopy with theoretical ab 
initio calculations. On the experimental side, imaging and electron-energy-loss spectroscopy 
(EELS) are combined via transmission electron microscopy (TEM) with sub-nm spatial 
and sub-eV energetical resolutions. On the theoretical side, the first-principles electronic- 
structure calculations are based on the local density functional theory (LDFT). This current 
progress of TEM and LDFT has been demonstrated for instance for the rhombohedral S7 
(1012) twin boundary in a-Al2O3,0Ji a model case of a realistic interface in a technologically 
relevant ceramic material. 

In the case of the rhombohedral twin boundary,^ both the quantitative analysis of high- 
resolution TEM (HRTEM) and EELS were able to discriminate between the lowest-energy 
interface model and another metastable one, which were predicted by ab-initio LDFT calcu- 
lations.! On the one side, for the high-resolution TEM combined with image simulation, the 
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discrimination was at the limits of the present experimental resolution and simulation ac- 
curacy. On the other side, the electron-energy-loss spectroscopy in a scanning transmission 
electron microscope (STEM) provided spatially resolved energy-loss near-edge structures 
(ELNES) in the oxygen i^-ionization edge, which, in combination with calculated local 
densities of unoccupied one-electron states, yielded a more clear distinction between the 
lowest-energy and metastable higher-energy interface models. 

In the present work, this effort in microscopic interface analysis is extended to the pris- 
matic S3 (1010) twin of a-Al203. In this case, the theoretical LDFT analysis results again 
in two metastable low-energy interface models. In distinction from the rhombohedral-twin, 
it is demonstrated that experimental EELS in a STEM with spatial resolution of 1—2 nm, as 
in Ref. ^ is not sufficient to discriminate the ELNES of the theoretically predicted models 
for the prismatic twin interface. It is predicted that for a new generation of analytical TEM 
instruments with spatial resolution of 0.1—0.2 nm, which are currently under development 
in academic and industrial TEM research laboratories,! the prismatic twin interface will pro- 
vide a promising benchmark for a successful achievement of atom-column-resolution ELNES, 
via a site-specific comparison of measured oxygen X-ionization edges and calculated local 
densities of states. A further challenge for quantitative high-resolution TEM image analysis 
is provided by the calculated interfacial arrangement of atoms. 

After a concise summary of the computational and experimental methods (Sec. II), we 
will report the search procedure for metastable variants of the twin boundary via empirical 
shell-model and ab-initio LDFT calculations (Sec. IIIA), the ab-initio LDFT analysis of the 
interfacial atomistics and energetics (Sec. IIIB), the calculated local densities of states and 
experimental ELNES (Sec. IIIC), and the simulated HRTEM images of the theoretically 
obtained two variants of the twin boundary (Sec. HID). The main conclusions will be 
summarized in Sec. IV. 



II. METHODS 
A. Computational approaches 

The theoretical analysis described in this work combines empirical atomistic modeling, 
based on an ionic shell-model potential J and first-principles electronic-structure calculations, 
based on the density-functional theorya in the local-density approximation for exchange and 
correlation! (local density functional theory, LDFT). 

The shell-model potential was parameterized for a;-Al203 by Lewis and Catlowi via 
fitting to experimental and theoretical data. The potential describes the crystal as fully 
ionic, with long-range electrostatic interactions (point charges on both atomic sites and 
charged shells on the anion sites), short-range overlap repulsion, and weak van der Waals 
interactions. The model reproduces quite well the structural properties of bulk and defective 
a-Al203, but fails to predict the energetic stability of the corundum structure relative to 
other possible polymorphs.!! In this work, the model is applied for a qualitative search of 
metastable translation states of a-Al203 (corundum) bicrystals, hence it is applied where 
its predictive power is well justified. The atomistic interface-structure variations with the 
shell-model potential are performed with the computer code GULP.t3 
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In our first-principles mixed-basis pseudopotential LDFT approach,llJllj the crystal va- 
lence wave functions are represented by a mixed basis set, consisting of plane waves and 
localized functions. In the specific case of AI2O3, the plane waves are limited by the maxi- 
mum kinetic energy Ep^=300 eV. Three localized functions with angular momentum / = 1, 
m = 0,±1, are centered at each oxygen site, and spatially confined by a sphere of radius 
Tio=l-lQ A. The interactions of the atomic nuclei and core electrons on the valence electrons 
are described with norm-conserving non-local ionic pseudopotentials. More details on the 
construction of the pseudopotentials, and on the predictions of the method for the structural 
parameters of bulk a-Al203 can be found in Ref. ^ The total energies for the prismatic 
twins are calculated by sampling the Brillouin zone of orthorhombic 120-atom supercells 
(see below) on a 2x1x1 k-point mesh. The forces on all atoms are explicitly calculated af- 
ter self-consistency of the electronic structure is achieved, and thus the method allows for 
efficient structural relaxation of the atomic positions. 

The crystal structure of a-Al203 (corundum or sapphire, space group R3c), projected on 
the (1210) plane, is shown in Fig. light and dark circles represent oxygen and aluminium 
ions, respectively, and the dashed line marks a prismatic (1010) plane. The twin boundaries 
are modeled with 120-atom orthorhombic supercells. These are defined by the lattice vectors 
ei = [1210], e2=[0001], and e3=[1010]. The supercell contains two grains of a-Al203 with a 
misorientation given by (1010) || (1010) and [1210] || [1210]: the three-dimensional periodic 
repetition generates a sequence of equivalent interfaces, separated by 6.8 A thick slabs of 
bulk corundum. 



B. Experimental procedures 

The prismatic twin bicrystal of a-Al2p3 was prepared by W. Kurtz, MPI fiir Metall- 
forschung Stuttgart, by diffusion bondingllB of two halves of highly pure single crystals of 
sapphire in ultra-high vacuum. The resulting bicrystal has a deviation of ?«0.1° from the 
perfect misorientation (measured by electron and X-ray diffraction). The impurity concen- 
tration at the interface was below the detection limit of energy-dispersive X-ray (EDX) spec- 
troscopy (0.3 atoms/nm^).tZl The sample preparation for transmission electron microscopy 
and spectroscopy involved conventional sawing, polishing, grinding, and finally ion milling 
at rather low energy of 2.5 keV, which resulted in highly radiation-damage resistant TEM 
samples.0 Electron energy-loss spectra (EELS) were recorded in a dedicated scanning trans- 
mission electron microscope (STEM, VG HB501UX) run at 100 kV and equipped with a 
parallel EELS spectrometer and an EDX detector. The measurements of spatially resolved 
energy-loss near-edge structures (ELNES) of the oxygen K edge were performed in a spot 
modus, with energy dispersion of 0.1 eV per channel and effective electron spot diameter 
between 1 and 2 nm. The energy resolution of the system was 0.8 eV. Extensive exper- 
imental TEM investigations of diffusion-bonded twin bicrystals of a-Al203 with various 



misorientations are reported in Ref. 17 



III. PRESENTATION AND DISCUSSION OF RESULTS 
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A. Metastable variants of the prismatic twin interface 

An initial search for metastable lateral translation states of the prismatic twin bicrystal 
was done by modeling the two grains as rigid blocks, and exploring the energy difference as 
a function of the relative rigid shift in the boundary plane (1010). A complete structural 
minimization was then performed only for the resulting equilibrium bicrystal configurations. 
The energetics was calculated via atomistic modeling based on an empirical shell-model 
potential.Q This simple but effective method has hitherto always succeeded in isolating the 
complete set of possible metastable translation states of the alumina bicrystals, ME which 
is the only goal of this part of the study. 

The lateral translation state T = Tiei + T2e2 is defined in terms of the lattice vec- 
tors ei = [1210] and 62 = [0001]. The reference bicrystal configuration, corresponding to 
(Ti, r2) = (0,0), was constructed from a block of bulk alumina by rotating one half of the 
crystal by 180° around [1010] with respect to the other half. Since the rotation axis contains 
an inversion point of the crystal, the rotation generates a mirror-symmetric twin boundary, 
denoted by M hereafter. Fig. |^ shows the projection of one grain on the prismatic (lOlO) 
plane: the solid lines correspond to the lattice vectors ei and 62 identifying the cell-edges of 
the orthorhombic supercell used in the calculations. The label M marks the position of the 
other grain in the reference M configuration, where the interfacial atoms of the two grains 
face each other at very close distance. 

The energetics of the bicrystal as a function of the translation state {Ti,T2) is shown 
in Fig. 0. Note how the highest energies correspond to the energy profile for Ti=0 (show- 
ing only two saddle points), while the lowest energies correspond to the energy profile for 
Ti=l/2. This can be easily understood by noting that the translation states with Ti=l/2 
reconstruct a highly ordered bulk-like atomic environment with hexagonal symmetry in the 
[0001] projection. Only two metastable boundary variants are identified, with translation 
states (Ti, r2) = (l/2,0) and (Ti, T2) = (l/2,l/3), denoted by G and S, respectively. In the 
most stable S variant, the two grains are related by a screw-rotation symmetry operation 
around e3=[10iO], while a glide-mirror symmetry operation with respect to (1010) relates 
the two grains in the metastable G boundary configuration. 

B. Energetics and atom arrangements 

The internal degrees of freedom of the two supercells describing the S and G boundaries 
were fully relaxed with respect to the LDFT atomic forces. The resulting interface structures 
are highly ordered and are shown in Fig. ^ and Fig. ^. The interface energy was calculated 
as the difference between the total energy of the relaxed supercell, and the total energy 
of an equal number of AI2O3 bulk formula-units, divided by the total interface area in the 
supercell (121.6 A^). 

The S boundary variant is the most stable structure of the prismatic twins, with an 
interface energy of 0.30 J/m^, which is the lowest interface energy of the alumina boundaries 
studied so far from first-principles (0.63 J/m^ for the E7 rhombohedral,! 0.73 J/m^ for the S3 
basal,@ and 1.88 for the S13 pyramidalEl twins). The metastable G variant is characterized 
by a higher interface energy, 0.49 J/m^, but still very low if compared to the basal and 
rhombohedral twins. These values correspond to zero excess volume (no expansion along T3), 
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therefore it is expected that a cell-volume minimization would further reduce the calculated 
interface energies. However, the energy lowering due to the inter-granular separation is 
always very small (of the order of some mJ/m^), particularly the one for highly ordered 
boundary structures. As a consequence, this is a case in which the relaxation with respect 
to the inter-granular separation will certainly be very small and can therefore safely be 
neglected. 

For completeness, the S and G boundary variants were also structurally relaxed according 
to the shell-model forces. The resulting structures were similar to the corresponding LDFT 
ones. The shell-model calculations predict the correct energy ordering between the S and G 
interfaces, with absolute values of 0.1 and 0.4 J/m^, respectively. It is worth noting that this 
is not a general rule, since the shell-model predicts the correct qualitative energy ordering 
between the variants of the rhombohedral S7 (Ref. ^ and of the S13 (Ref. ^) twins, but 
fails in the case of the basal S3 twin (Ref. p!8[) . 

1. The screw-rotation twin S 

The structure of the prismatic screw-rotation twin interface relaxed via LDFT is shown 
in Fig. ^. The interface, marked by dashed lines, is decorated by Al and O atoms. It is 
highly ordered, so much that the projection of the twin on the (0001) plane (upper panel) 
is indistinguishable from the one of the corundum single crystal (upper panel of Fig. |l|). If 
the small interfacial relaxation of the oxygen sublattice is neglected, the boundary does not 
affect the symmetry of the anion sublattice, that is hexagonal close-packed throughout the 
supercell. In fact, the twin interface arises mainly from a modified packing of the sublattice 
of cations, which fills two thirds of the octahedral interstitial sites in the anion sublattice, 
leaving one third of the sites empty (denoted as "voids" in the following). The arrangement 
of the voids across the twin boundary can be clearly identified by the symmetric V-shaped 
pattern. 

The Al-0 bond lengths of the interfacial atoms are included in Table |, together with 
the Al-0 distances in bulk q;-A1203 (short bs=1.84 A and long bi=1.96 A bonds). The 
interfacial All and A12 atoms are six-fold coordinated, with three bond lengths centered 
around b^ and three around b;, hence the cation sites are very bulk-hke. Since the boundary 
affects mainly the packing of the cation sublattice in the boundary region, the main changes 
in coordination number and local atomic environment are shown by the interfacial oxygen 
atoms. These are distributed between an equal number of three-fold under-coordinated 
(03), four- fold coordinated (04), and five-fold over-coordinated (05) atomic sites. The 04 
anions are bulk-like tetrahedrally coordinated, with two short bonds centered at hg and two 
long bonds centered at b/. The under-coordinated 03 oxygens have three short bonds in a 
planar geometry, while the 05 anions have one short and four long bonds. 

The atoms in the first atomic layer away from the interface are already bulk-like, with 
Al-0 distances of 1.84 and 1.96 A, therefore showing that the perturbation in the corundum 
crystal symmetry due to the presence of the boundary is extremely confined to the interface 
plane, with a transverse extension of less than 2.5 A (compared for instance to 5-6 A in the 
E13(1014) twinS). 
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2. The glide-mirror twin G 



The structure of the prismatic ghde-mirror twin relaxed via LDFT is shown in Fig. |^. 
The dashed hne marks the position of the interface, that bisects two planes of atoms in 
this case. The boundary structure has many similarities with that one of the S boundary: 
the (0001) projection is bulk-like and the interface is determined by the cation sublattice 
only. As a result, the boundary cations are six- fold bulk-like coordinated, while the main 
changes in the local environment are shown by the boundary anions. Similarly to the S 
variant, these are distributed between three-fold under-coordinated (03), four-fold bulk-like 
coordinated (04), and five- fold over-coordinated (05) sites (Table ^. 03 has three short 
bond lengths, 04 has two short and two long ones, and 05 has one short and four long ones. 
The inner atoms neighboring the boundary layer are in a bulk-like environment, therefore 
the transverse extension of the boundary is again very small in this case (< 4 A). 

The main distinctions between the two twin models can be summarized by these two 
points: (1) The stacking sequence of the inequivalent boundary anions along the [0001] (03- 
04-05 in the S boundary, 05-04-03 in the G boundary), which is accessible by ELNES 
measurements with spatial single-atom-column resolution. (2) The pattern formed by the 
cation-void sites across the boundary (V-shaped with one single void at the S boundary, V- 
shaped with a pair of voids at the G boundary), which is accessible by quantitative HRTEM 
analysis. 

3. The free prismatic surface 

The free prismatic surface was relaxed by means of LDFT in a 60-atom supercell con- 
taining six layers of bulk alumina. The transverse size of the supercell was chosen to be big 
enough so that to minimize the spurious interaction between the free surfaces. These were 
separated by 9.5 A of vacuum and by 6.8 A of bulk corundum. 

The relaxation involves mainly the two outermost atomic layers, with minor atomic 
displacements (< 0.05 A) in the inner third layer. The main features of the relaxed surface 
are the outward expansion (0.2 A) of the external anions, the inward contraction (0.3 A) of 
the external cations, and the large lateral displacement (0.4 A) of the surface anions. The 
calculated LDFT value for the surface energy is 2.58 J/m^. 

C. Projected densities of states and ELNES of oxygen 

Experimentally, the boundary structure of the E3 prismatic twin was studied by EELS, 
measuring the K-edge ELNES of oxygen. When ELNES is analyzed by means of the spatial 
difference technique,Ey it gives an indirect information about the local bonding and coordina- 
tion through a direct measure of the unoccupied electron states present in the neighborhood 
of an excited atom. These electron states can be obtained from the calculation of the double 
differential scattering cross section for the excitation of an atom by inelastic electron scatter- 
ing.i3 In the special case of the oxygen Is — *• 2p excitations, and in the dipole approximation 
for the electron scattering, the scattering cross section, and therefore the ELNES K-edge 
for an oxygen site, can directly be related to the calculated site- and p-projected density of 
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states (p-PDOS, p denotes the angular momentum /=1).E30 Final-state (core-hole) effects 
are included in these calculations via the Z+1 approximation^where the excited oxygen 
atom is simulated by a substitutional fluorine impurity atom.E3c3 The p-PDOS at different 
excited oxygen sites are calculated from the LDFT one-electron eigenvalues on the 2x1x1 
k-point mesh, and smoothened by a Gaussian of 0.8 eV to mimic instrumental and life-time 
broadening effects. 

To see whether the calculated lowest-energy boundary structure is the same as the ex- 
perimentally observed one, we compare the calculated p-PDOS of the S and G interfaces 
with the measured ELNES K-edge in Fig. ^. The experimental interface spectrum shows a 
pronounced peak at 541 eV with a shoulder at higher energies, and a minor peak around 560 
eV. The calculated p-PDOS curves for the S and G interfaces are shifted along the energy 
axis to match the first peaks above the edge onset with the peak at 541 eV in the experimen- 
tal spectrum. The calculated p-PDOS are very similar to the experimental ELNES, with 
a good agreement between the relative position and relative height of the peaks. However, 
they are also very similar to each other, so that, from this comparison, it is not possible 
with sufficient confidence to discriminate one or the other theoretical boundary structure 
from the real interface. 

The similarity between the jo-PDOS can be understood by recalling that they reflect the 
S and G boundary atomic environments, which are qualitatively and quantitatively very 
similar (Table |). Fig. |^ shows the breakdown of the total interfacial jo-PDOS, displayed 
in Fig. 1^, into the individual contributions corresponding to the 03, 04, and 05 anions. 
The difference in coordination number affects mainly the low energy region of the spectrum. 
Over- and under-coordination drives the distinctive "peak and shoulder" shape of the bulk- 
like atoms, towards a single peaked, and double peaked shape, respectively. These distinctive 
signals coming from three-, four-, and five-fold coordinated atoms should in principle be 
measurable by EELS in a STEM with spatial single-atom-column resolution. In this case it 
would indeed be possible to discriminate between the two boundary structures, even though 
their p-PDOS are nearly indistinguishable. A line scan along the boundary layers would 
produce periodic sequences of spectra corresponding to the 03, 04, and 05 anions. The 
atom-column specific signals for the S and G variants would be similar, but their spatial 
sequence along the direction [0001] parallel to the interface would be reversed in order, thus 
allowing to discriminate between the two interfaces: ... -03-04-05-... in the case of the S 
boundary, ... -05-04-03-... in the case of the G boundary in bicrystals with equivalent bulk 
orientations. However, for the time being such a level of spatial atom-column resolution is 
set as a challenge for experimental STEM-EELS microanalysis.! 

D. Simulated HRTEM images 

For a bicrystal sample of q;-A1203 illuminated by a spatially extended, parallel electron 
beam along the [12l0] zone axis in a high-voltage TEM with atomic-resolution capability, 
the main feature of the resulting HRTEM image is typically a strong intensity-contrast 
pattern with intensity extrema at the positions of unoccupied cation sites (cation void sites) 
in the corundum structure of a;-Al203. The detailed atom-column resolved projection of 
the bicrystal structure is not directly accessible from the measured HRTEM micrograph. It 
can be deduced by comparing it via quantitative matching with simulated HRTEM images 
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from atomistic- model candidate structures-Ejo For optimum agreement, all atom-column 
positions as well as imaging parameters (e.g., defocus, specimen thickness, lens aberrations 
and transfer function of the recording medium) need to be determined variationally. 

As discussed above, we have extracted two metastable structures which should in prin- 
ciple result in different HRTEM images in the (1210) projection. The simulated HRTEM 
images of the S and G supercells corresponding to a TEM sample thickness of 6 nm and 
to a defocus value of —20 nm, which are typical parameters for HRTEM of a-Al^Os in an 
instrument like the JEM ARM-1250 at the MPI fiir Metallforschung Stuttgart ,iHil are dis- 
played in Fig. 1^. The two images show characteristic features that allow a discrimination 
between the two interfaces. In the S variant, the white spots located at the boundary plane 
are visually equivalent to the ones inside the bulk grains. In the G variant, the interfacial 
white spots have a distinctively oversized spot, arising from two neighboring spots merging 
together. The difference in the micrographs does not arise from the interfacial translation 
states, since both the S and G variants produce a symmetric V-shaped pattern of octahedral 
void sites, but from the local boundary structure. In particular, the rhombohedral planes of 
voids in different grains intersect the interface plane precisely at one cation-void site in the 
S variant (Fig. and on an oxygen basal plane in the G variant (Fig. This is schemat- 
ically shown by the dashed lines in Fig. ^, where the cation-void sites are associated to the 
white spots. These images suggest that an high-resolution micrograph would discriminate 
between the two boundary variants. Experimental HRTEM imagesSS show a contrast 
pattern compatible with the simulated one for the S interface, apparently confirming the 
LDFT prediction based on the lowest interface energy. Hence, for such experimental im- 
ages a quantitative HRTEM analysiscHta starting from our calculated two interface models 
appears promising (cf. Ref. p. 



IV. CONCLUSIONS 

In this work, the energetic stabilities, atomistic arrangements and electronic states at the 
prismatic S3 (1010) twin interface in q;-A1203 were analyzed theoretically by atomistic shell- 
model and ab-initio LDFT calculations, and experimentally by spatially resolved STEM- 
EELS measurements of the oxygen i^'-ionization-edge ELNES. 

Two metastable microscopic variants for the prismatic twin with low interfacial ener- 
gies resulted from the theoretical analysis. They exhibit subtle differences in the local ar- 
rangements of atoms at their interfaces, which should be distinguishable via TEM bicrystal 
experiments. 

Experimentally, it was demonstrated that a state-of-the-art STEM instrument, with 
spatial resolution of 1—2 nm and energetical resolution of about 0.8 eV, is not yet capable 
to resolve the subtle differences in the interfacial ELNES of this very bulk-like twin interface 
of a-Al203. However, with the next generation of analytical TEM instruments, which are 
designed to achieve resolutions of 0.1—0.2 nm in space and about 0.2 eV in energy and 
currently approaching the stage of operation, the two theoretical twin variants will become 
distinguishable. 

Hence, it is proposed that the present theoretical ab-initio results for the prismatic twin in 
a;-Al203 provide a very promising benchmark case for the coming analytic TEM instruments: 
(a) to see whether single-atom-column resolution is really achieved by detecting significant 
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differences in tlie oxygen X-cdgc ELNES for tlie tliree distinct interfacial oxygen sites; (b) to 
clieck wlicthcr the ab-initio LDFT prediction of the S variant with lowest interface energy is 
also the experimentally observed case. For the ELNES experiment, the spatial sequence of 
the three oxygen sites is decisive. Alternatively, for a HRTEM experiment combined with a 
quantitative image-matching analysis, the contrast patterns originating from arrangements 
of cation- void sites across the interfaces discriminate the two theoretical variants. 
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TABLE I. Equilibrium Al-0 bond lengths (in A) in the inter facial regions of the S and G 
supercells. For comparison, the Al-0 distances in bulk Q;-Al203 are also included. Atoms are 
numbered as in Fig.s 4 and 5 
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FIG. 1. Corundum structure of Q;-Al203 viewed from the [1210] direction. Light-gray and black 
circles represent oxygen and aluminium atoms, respectively. The dashed line marks the (1010) 
plane leading to the prismatic S3 interface. 
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FIG. 2. Atomic structure of one grain projected on the prismatic plane (1010). The in-plane 
extension of the computational cell is marked by the lattice vectors ei = [1210] and e2 = [0001], 
and by solid lines. M, G and S denote the translation states of the upper grain, leading respectively 
to Mirror, Glide-mirror and Screw-rotation twins. The atoms are represented as in Fig. [l|. 
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FIG. 3. Shell-model total energy (Ry/supercell) for the twin interfaces as a function of the lateral 
translation state Tiei + T2e2, with ei = [i2i0] and 62= [0001]. 
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FIG. 4. LDFT result for the relaxed supercell of the S interface viewed along the [0001] (above) 
and the [1210] (below) directions. The supercell contains two equivalent interfaces (dashed lines), 
one in the middle, the other at the extremes arising from the periodic boundary conditions in the 
[1010] direction. The atoms are represented as in Fig. [l|. 
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FIG. 5. LDFT result for the relaxed atomic structure of the G interface viewed along [0001] 
(above) and the [1210] (below) directions. Atoms, labels, and symbols as in Fig. 0. 
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FIG. 6. Comparison between the measured ELNES oxygen JC-edge signal and the calculated 
p-PDOS for the interfacial anions in the G and S boundary variants. (The p-PDOS is the sum of 
the individual p-PDOS for the three interfacial oxygen sites 03, 04, and 05). 
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FIG. 7. Individual contributions D(E) of the anions 03, 04, and 05 to the total interfacial 
oxygen p-PDOS calculated for the S (top) and G (bottom) boundaries. 
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